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<^ ■ Abstract 

, For construction and classification of tlie natural integrable systems we propose 

r— 1| to use a criterion of separability in Darboux-Nijenhuis coordinates, which can be 

ly^ • tested without an a priori explicit knowledge of these coordinates. 

d ■ 
•1—1 ' 

'a ■ 1 Introduction 

CN I The search of integrable dynamical systems is one of the most fascinating branch of 

^ ■ classical physics. Integrable systems are quite rare and still only a few examples are 

■ known. In this paper we will present a new direct method for construction of natural 
^ , integrable systems on tuN bi-hamiltonian manifolds. 

I A bi-hamiltonian manifold is a smooth manifold M endowed with a pair of compatible 

^ I Poisson tensors P and P' associated with the Poisson brackets {., .} and {.,.}' respectively. 

■ The class of manifolds we will consider are particular bi-Hamiltonian manifolds where one 
• S ■ of the two Poisson tensors in non degenerate and thus defines a symplectic form tu = P~^ 

^ , and a recursion operator = P'P^^. 

^ \ Dynamical integrable systems on M with functionally independent integrals of motion 

Hi, ... , Hn in the bi-involution 

H ■ 

.?^.\ {H,,Hj} = 0, and i/^}' = 0, for alH,i (1.1) 

will be called bi-hamiltonian systems for brevity. 

A Hamiltonian system is natural if its Hamiltonian is a sum of a positive-definite 
kinetic energy and a potential. For instance, on the symplectic manifold M ~ M^", i.e. 
in the Euclidean space, the natural Hamiltonian is equal to 



H = J2p^ + V{q). (1.2) 



i=l 



In the known direct method to find natural integrable systems JU] we substitute natural 
Hamiltonian and some ansatz for other integrals of motion into the equations (jl.lj) and 
try to solve these equations. 

According to [HI EI integrals of motion Hi, ... , Hn are in bi-involution (jl.lj) if they 
are separable in Darboux-Nijenhuis coordinates which are canonical with respect to to 
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and diagonalize A^. So, instead of classification of bi-hamiltonian systems we can classify 
Hamiltonians whose associated Hamilton- Jacobi equations can be solved by separation of 
variables in Darboux-Nijenhuis coordinates. 

A test for separability of a given Hamiltonian H in a given coordinate system was 
proposed by Livi-Civita in 1904 jTn|. Using Darboux-Nijenhuis coordinates associated 
with the Livi-Civita criteria may be rewritten in the following form p[] 

d{N*dH){VH,VH) = 0. (1.3) 

Here Vh is distribution generated iteratively by the action of on the vector field X^- 
The relevance of a new form of the Livi-Civita criteria is that condition p.3p can 
be checked in any coordinate system. So, for a given Hamiltonian H we can look at 
fjl.3|) as one of the equations to determine torsion free tensor A^ and hence the separated 
variables. This scheme, in its basic features, has already been considered in the literature 
and applied to various systems, see [H El 13 1121 UZl UHl UHl jSHj and references within. 

In this paper we will take a different logical standpoint. We propose to solve equation 
fll.3|) with respect to a torsionless tensor A^ and a natural Hamiltonian H simultaneously. 
We will prove that it is indeed possible and, actually, easy to solve such an equation 
by means of a couple of natural Ansatze. As an example, we reproduce a lot of known 
natural integrable systems on the plane with cubic and quartic integrals of motion pUj- 
The corresponding Poisson pencils and tensors A^ are new. 

2 Bihamiltonian geometry and separation of variables 

In this section we sketch the main points of the separation of variables theory for bi- 
hamiltonian manifolds, referring to d for complete proofs and more detailed 
discussions. 

2.1 uN manifolds 

The phase spaces of integrable systems will be identified with cjA^-manifolds IZl or 
Poisson-Nijenhuis manifolds ^lEl endowed with a symplectic form u and a torsion free 
tensor A^ satisfying certain compatibility conditions. 

Definition 1 An uN manifold is a bi-hamiltonian manifold in which one of the Poisson 
tensor (say, P) is non degenerate, i.e. det P ^ . 

Therefore, M is endowed with a symplectic form u = P^^ and with a tensor field 

N = P'p-\ (2.1) 

which is recursion operator (formally). The operator A^ is called a Nijenhuis tensor J7j 
or a hereditary operator |H| as well. 

In order to construct operator A^ by the rule ()2.H) we have to get a pair of compatible 
Poisson tensors P and P'. 

Definition 2 The Poisson tensors P and P' are compatible if every linear combination 
of them is still a Poisson tensor. The corresponding linear combination P^ = P + AP', 
X E M., is a Poisson pencil. 
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Recall a bivector P' on M is said to be a Poisson tensor if the Poisson bracket 
associated with it 



2n 



{f{z),g{z)y=J2P:ki 



df{z) dg{z) 



., , dzi dzk 

satisfies the Jacobi identity, that is, if for all f,g and h G C°^(M) 

{/, {g. h]'}' + {g, {K fYY + {K {/, g}'}' = o . (2.2) 

Here z = {zi, . . . , Z2n) is a point of M in some coordinate system. 

From the Jacobi identity we can get a system of equations on the entries P/kiz) of 
the tensor P'. So, in order to construct the Poisson tensor P' we have to solve algebraic 
equations 

P^j 'Pji'' ^'^ 1; • • • ) 2?T., 

and partial differential equations 

^" / dP' BP' dP'\ 

Y.\Prm^ + P;^^ + PL^]=^, 2,j,fc = l,...,2n, (2.3) 

yj^—-^ \ ^Zfn (yZjYi OZfji J 



with respect to some unknown functions Plj^{z). 

In order to construct the Poisson pencil on uN manifold associated with a given 
tensor P we have to solve equations ()2.3j) simultaneously with the equations 

/ dP^ dP^ dP^\ 

Y.{^'^^ + P'rngf + PLQf]=0, ^, J, /c = 1 , . . . , 2n, (2.4) 



m=l 



which follows from the Jacobi identity for the Poisson bracket {., .}a = •} + .}', 
where A is an arbitrary numerical parameter. 

For any solution P' of the equations ()2.3II2.4|1 there is free torsion tensor 1)2.11] . 

Theorem 1 The Nijenhuis torsion of N, 

Tn{X, Y) = [NX, NY] - n(^[NX, Y] + [X, NY] - N[X, Y]^ = 0, (2.5) 

vanishes as a consequence of the compatibility between P and P'. 

In a given coordinate system on M entries of the Nijenhuis torsion of are equal to 

^ fdNi dm dNY" ■ dNr ■ 

TJN) = V — -^iV™ iNr + — —N' - —J^N' 



m=l 



In general there are infinitely many solutions of the system of equations ()2.3H2.4p or ()2.5|) . 

In order to identify 2n separated variables with Darboux-Nijenhuis coordinates we 
have to consider a special class of uN manifolds. 

Definition 3 A 2n- dimensional uN manifold M is said to be semisimple if its recursion 
operator N has, at every point, n distinct eigenvalues Ai, . . . , A„. It is called regular if 
eigenvalues of N are functionally independent on M . 

Of course, equations ()2.3II2.4|) also have infinitely many solutions on semisimple regular 
uoN manifolds. In order to get a bounded set of the solutions we have to fix form of the 
admissible Poisson tensors P' . 
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Example 1 Let us consider canonical Poisson tensor in the Euclidean space 

^= ( ly I = diag(l,...,l) 
associated with the usual Poisson bracket for arbitrary functions f,g& C°°(M) 

\dpi dqi dqi dpi J 

At n = 2 there are two nontrivial solutions of the system of equations ()2.3II2.4|) 
P' 



/o 












/o 





ai - qf 


-qiq2 \ 


* 





Qi 


2q2 


and P' = 


* 





-qiq2 


«2 - ql 


* 


* 





-Pi 




* 


* 





-qiP2 + q2Pi 




* 


* 


J 




V* 


* 


* 


/ 



(2.6) 



(2.7) 



(2i 



which are matrix polynomials in variables z = {p, q). The Nijenhuis torsion of the corre- 
sponding recursion operator fl2.1|) vanishes and its eigenvalues are distinct in the both 
cases. 



2.2 The Darboux— Nijenhuis coordinates 

According to (3j the recursion operator A^ exists for overwhelming majority of integrable by 
Liouville systems. Therefore classification of the Poisson pencils has to coincide with the 
classification of integrable systems, which has a more rich history [lOj. This connection is 
based on the separation of variables method in which separated variables for the Hamilton- 
Jacobi equation are identified with Darboux-Nijenhuis coordinates |16j . 

Definition 4 A set of local coordinates {xi,yi) on an uN manifold is a set of Darboux- 
Nijenhuis coordinates if they are canonical with respect to the symplectic form 



u = P ^ = dyi A dxi , 

i=l 

and put the recursion operator N in diagonal form 

^ = 5Z [d^. ®dxi + — dyi J , (2.9) 

This means that the only nonzero Poisson brackets are {xj,?/^} = 5ij and = Aj5.y. 

According to ()2.9|) differentials of the Darboux-Nijenhuis coordinates span an eigenspace 
of the adjoint operator N* 

N*dxi = Xidxi, N*dyi = Xidyi . (2.10) 

As a consequence of the vanishing of the Nijenhuis torsion of A^, the eigenvalues Aj always 
satisfy 

N*d\i = XidXi . 

It allows us to find very easy one special family of Darboux-Nijenhuis coordinates. 
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Theorem 2 /n a neighborhood of any point of a regular uN manifold where the 
eigenvalues of N are distinct it is possible to find by quadratures n functions /ii, . . . 
that, along with the eigenvalues Ai, . . . , A„, are Darboux-Nijenhuis coordinates. 

In this case the coordinates Xj can be computed algebraically as the roots of the minimal 
polynomial of 

A(A) = (det(iV- AI)) =JJ(A-Aj). 

On the contrary, the conjugated momenta must be computed (in general) by a method 
involving quadratures [HI E] • 



Example 2 For the first Poisson tensor in ()2.8|1 the minimal characteristic polynomial 
of is equal to 

A(A) = A^ - g^A - I = (A - Ai)(A - As) . 
The corresponding special Darboux-Nijenhuis coordinates 



Ai,2 = ^\^q2±\/ql + qlj , /ii,2 = P2 - — \q2 ^ \j ql + ql 

coincide with the parabolic coordinates on the plane and their conjugated momenta. 
For the second Poisson tensor in ()2.8|1 the corresponding minimal polynomial 



A(A) = A^ + {ql + ql - ai - 02) A + aia2 - a2ql - aiqj = (A - Ai)(A - A2 
is a generating function of the elliptic coordinates Ai_2 on the plane |2]. 



2.3 Integrable systems in the Jacobi method 

The Darboux-Nijenhuis coordinates (xj, yi) or (Aj, /Xj) may be identified with the separated 
variables for a huge family of integrable systems. In order to construct these integrable 
systems in framework of the Jacobi method we can use the following theorem. 

Theorem 3 Let [Xi, are canonical coordinates {Xi, fij} = 6ij. The product of n one- 
dimensional Lagrangian submanifolds 

Cj : $j(Aj, /ij, «!, ...,«„) = with det 

is an n-dimensional Lagrangian submanifold jF(a) = Ci x ■ ■ ■ x C„,. The solutions of the 
separated equations \2. 1 1]) with respect to n parameters au are functionally independent 
integrals of motion Hk = ak{X, fi) in involution. 

The main problem in the Jacobi approach is to choose separated equations $j = for 
which the Hamilton function H = Hi has a natural form ()1.2|) in the initial physical 
variables (p, q). 

We can reformulate this problem by using the separability condition ()1.3|) which im- 
plies that the distribution is integrable. So, there exist n independent local functions 
i^i, . . . , Hn that are constant on the leaves of 'Dh- 



7^0 



(2.11) 
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Theorem 4 |^[^Let M be a semisimple regular uN manifold. The functions Hi, ... , Hn 
on M are separable in Darboux-Nijenhuis coordinates if and only if they are in bi- 
involution Al.l]) . The distribution D tangent to the foliation defined by Hi, . . . , Hn is 
Lagrangian with respect to u and invariant with respect to N. 

The invariance of the Lagrangian distribution T> means that there exists a matrix F with 
eigenvalues (Ai, . . . , An) such that 

n 

N*dHj = Y,FjkdHk, i = l,...,n. (2.12) 

fc=i 

Here F is the control matrix with respect to integrals of motion Hi, ... , Hn, which form 
the Nijenhuis chain jHllZl- 

Using this theorem we can formulate a new method to find natural integrable systems. 
Namely, substituting a torsionless tensor A^, some control matrix F and a natural Hamilto- 
nian H = Hi fll.2|) into the n relations ()2.12|) . one gets a system of n differential equations 
on the potential V{q) and on the remaining integrals of motion H2, ■ ■ ■ , Hn. Solving these 
equations one gets a natural integrable system separable in Darboux-Nijenhuis coordi- 
nates. In order to avoid a lot of intermediate calculations we can solve equations ()2.12j) 
together with equations ()2.3|2.4|) . 

The proposed construction inherits the main disadvantage of the Jacobi method be- 
cause we have to choose control matrix F for which equations (|2.12j) have non-trivial 
solutions. 

Example 3 If Hk = j^tiN^ then the control matrix in ()2.12|) is equal to 

/O 1 ■■■ o\ 
1 ■■■ 



; ■■■0 1 

CTn+l ■■■ Oij 



(2.13) 



Here (7^ are the elementary symmetric polynomials of degree k on the eigenvalues of A^. 
In this case relations (|2.12p are the Lenard-Magri recurrence relations 

N*dHi = dHi+i or P'dHi = PdHi+i, (2.14) 

such that the corresponding vector fields Xh^ are bi-hamiltonian jT^l EI • 

Remark 1 It is known that for finite-dimensional systems the bi-hamiltonian property of the 
vector fields is a very strong condition and too restrictive for construction of natural integrable 
systems jTl E]. Of course, we could easy solve equations 1)2. 12() with such matrix F 1)2. 13() . 
but for the majority of the known Poisson pencils system of equations (|2.12j) or (|2.14|) is an 
inconsistent system if Hi = H is a natural Hamiltonian ()1.2() . 



2.4 The uniform Stackel systems 

The equations ()2.12j) do not give explicit information on the form of the separated equa- 
tions (j2.1H) . In order to restrict a set of the possible separated equations $j = we will 
consider the following family of the separated equations 

n 

^i{fj,i,\i,ai, . . . ,an) = ^Sij{\i) aj - Si{^i,Xi) = , Hi = ai , (2.15) 

3=1 
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where is a non-degenerate Stackel matrix and s is a Stackel vector. The entries Sij and 
Si depend on a pair of the separated variables Aj and /ij only. 

Theorem 5 |3 Let Hi, . . . ,Hn are independent integrals of motion in bi-involution 
m.l]) defining a bi-Lagrangian foliation on a regular semisimple ujN manifold. If the 
control matrix F satisfies 



N*dF = FdF, that is N*dFij = ^ FikdFt 



kji 



i,j = l,...,n, (2.16) 



k=l 



then the functions Hi, ... , Hn are Stackel separable in the Darboux-Nijenhuis coordinates. 
In this case the Stackel matrix 5" (j2.15|) is defined by 

F = S~^AS, A = diag(Ai,...,AO. 

Construction of the corresponding Stackel vector s is not so algorithmic. 

This theorem could be very useful for a construction of natural integrable systems 
of the Stackel type. Namely, substituting a (1,1) tensor and a natural Hamiltonian 
H = Hi ()1.2|) into the relations ()2.5|) . ()2.12|) and ()2.16|) one gets the closed system of 
algebro-differential equations 



T, 



N 



0, 



N*dH = FdH, N*dF = FdF, 



(2.17) 



on potential V{q), integrals of motion H2{p, q), . . . , H^ip, q), entries Nij{p, q) of the recur- 
sion operator and entries Fij {p, q) of the control matrix. 

Unfortunately, we can not solve equations ()2.17|) in general. More precisely, usually 
equations ()2.17|) have infinitely many solutions. The main reason is that we fix form of 
the Hamiltonian H ()1.2|) . but an admissible set of the Stackel vectors s in ()2.15|) remains 
undefined and, therefore, unbounded set. So, in order to solve equations (|2.17p in practice 
we have to narrow a class of the separated equations (|2.15|) once more. 

Below we will determine a class of the separated equations by using the Stackel matrix. 
Recall the Stackel matrix is a n x n block of the transpose Brill- Noether matrix, which is a 
differential of the Abel- Jacobi map associated with a product J-'{a) of the algebraic curves 
Ci (j2.1ip . see 1^ and references within. There is one to one correspondence between j-th 
rows of the Stackel matrix S and a basis of the differentials on the j-th curve Cj. 

Let us consider uniform Stackel systems j21J| for which the Lagrangian submanifold 
jF(a) is the n-th symmetric product of a hyperelliptic curve C. For the standard basis of 
the holomorphic differentials on C associated with the Veronese map the corresponding 
Stackel matrix is equal to 



S 



Ai l\ 



V„ 1/ 



n-l 



(2.18) 



It is easy to prove that the control matrix 



F = S~^AS 



(ai I 
a2 1 



cr„ ■ • • 



fe=0 



... 

1 
0/ 



(2.19) 
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fulfills inilZI- Here aj are elementary symmetric polynomials on eigenvalues Xj of 

tensor 

A(A) = A" - (TiA"-^ - aaA"-' a, . (2.20) 

Solving the separated equations ()2.15|) with respect to ai, . . . ,a„ one gets the integrals 
of motion Hi as functions on the separated variables 



n 



The corresponding vector fields Xh- are tangent to a bi-Lagrangian foliation, but they 
are not, in general, bi-Hamiltonian. These vector fields are Pfaffian quasi-bi-hamiltonian 
fields in the terminology of 



Remark 2 The Stackel matrix S H2.18() is one of the most studied matrices, which appears 
very often in various applications j2i |6| |?| 1121 121j . But we have to keep in mind that there are 
many other Stackel matrices associated with hyperelliptic or non-hyperelliptic curves. 

Example 4 Let us consider natural integrable systems on the plane separable in the 
polar coordinates (r, 4>,Pr,P<j>)- The corresponding second Poisson tensor P' is degenerate 



(0 





/p'2 


o\ 


* 











* 


* 








V 


* 


* 





(2.22) 



Operator ()2.H) has two different eigenvalues Ai = and A2 = 0. Integrals of motion 



Hi=Pr + ^ + Vir), H2=p^ (2.23) 



give rise to the following control matrix F and the separated equations $1,2 








P0 - ^2 = 



The corresponding Stackel matrix S is non-homogeneous in contrast with ()2.18|) . If V{r) is 
polynomial then the separated equations $1 = determines hyperelliptic curve, whereas 
curve associated with $2 = is non-hyperelliptic. 

Summing up, the separation of variables theory for uoN manifolds outlined above 
provides intrinsic and algorithmic recipe to construction of the natural integrable systems 
of the Stackel type. 

The algorithm consists of solution of the extended system of equations ()2.3II2.4|) and 
()2.12p with some fixed Stackel matrix 5*. These equations contain the following unknown 
functions Pj^-, V and H2, . . . , Hn- Solving this overdetermined system of equations in phys- 
ical variables we get a Poisson pencil, recursion operator, Darboux-Nijenhuis coordinates 
and the corresponding natural integrable system simultaneously. 

Remark 3 In the proposed algorithm we have to substitute some concrete Stackel matrix S 
into the equations (|2.12() and then to solve these equations. In general matrix S defines only first 
part of the separated equations (|2.15|) whereas the Stackel vector s remains arbitrary. However 
if we use the Stackel matrices ()2.18() associated with the Abel-Jacoby map on the hyperelliptic 
curves, we implicitly fix the Stackel vector s as well, because the separated equation has to 
generate namely hyperelliptic curve. 
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3 Bi-hamiltonian systems on the plane 



In this section we substitute natural Hamiltonian H ()1.2|) and some polynomial anzats for 
the Poisson tensor P' in (|2 .312 .412 .121) and describe nontrivial physical solution of these 
equations by ri = 2. 



Definition 5 Solution of the equations \2/J\2.4\^2.1^/3^) is nontrivial, if the functions Hi 
are functionally independent integrals of motion, dHi ^ for any i, and potential V{q) is 
some real function, V{q) ^const. 

Moreover, throughout the rest of the section we consider only physical solutions for which 
all the integrals of motion are polynomials in momenta. 

Example 5 Substituting an independent on momenta tensor P' and the following Hamil- 
tonian Hi = H = P1P2 + V{qi, 52) into (j2.3|2.4f2.12|l we can get the following solution 



P' 



/O 0\ 

* fiqi) 

* * 

* * oy 



Hi=pm + i^,+f{qi), 

H2 = Jqi exp{Cipl) . 



It is example of the non physical solutions which will not considered below. 

Of course, integrable systems with two-degrees of freedom are one of the well studied 
integrable systems and, therefore, we can only obtain some new tensors N for known 
integrable systems listed in ^2] • Nevertheless below we present two new integrable systems 
on the plane with pseudo-Euclidean metric as well. 



3.1 Linear Poisson tensors and the L-systems 



Let Q be an n-dimensional Riemannian manifold endowed with a conformal Killing tensor 
L of gradient type with vanishing Nijenhuis torsion According to ^21 its cotangent 
bundle M = T*(5 is an ujN manifold, whose recursion operator 



'dq' 



Qqk Qq'i j Qp 



— , N— = yL'— 

' dpk ~^ * dpi 



i=l ^ I] 

is a complete lifting of operator L. In this case the Poisson tensor 



(3.1) 



P' = NP 



( ° 

Li 

\ 




(3.2) 



is a linear matrix polynomial in momenta. 

Substituting a Poisson tensor P' fl3.2|) and a natural Hamiltonian Hi = T + V ()1.2p 
into the system of equations ()2.12|) and expanding it in powers of the momenta we can 
see that ()2.17j) splits into the equations on L and on V 



diCxrd' -Tdai) =0, d{Cxv 0' - Vdai) = 
and the cyclic recurrence relations for other integrals of motion 

A^* dHi = dHi+i + OidHi , i = 2,...,n, 



0, 



(3.3) 



(3.4) 
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which we have to solve starting with i = n and i^„+i = 0, see e.g. [U [T21 12S1- 

Here £ is a Lie derivative along the vector field Xt or Xy and 0' = ^"^=1 U^pidq^ is 
a deformation of the Liouville form = ^Pjdq^ for any set of fibered coordinates (p, q). 

In |23| 125 the computer program for finding of the L-tensors and the associated 
integrable systems was constructed. This software allows to solve equations (|3.3|) and 
(j3.4|) on the Riemannian manifolds with a positive-definite metric. 

The eigenvalues on ()3.1|) coincide with the eigenvalues Ai(q'i, . . . , g„) of the operator 
L defined on the configuration space Q. As a consequence, functions ai depend on coordi- 
nates q only and, according to ()3.4|) . for the natural integrable systems with Hi = T + V 
()1.2|1 all the integrals of motion H2, . . . Hn are quadratic polynomials in the momenta 

The associated with 1)3.11) Darboux-Nijenhuis coordinates are related with initial 
physical variables [p, q) by the point canonical transformations. Indeed in this case 
Darboux-Nijenhuis coordinates coincide with well-studied orthogonal curvilinear coor- 
dinates on the Riemannian manifolds. 

Example 6 On the plane Q = ffi^ there are four L-tensors, which give rise to the natural 
integrable systems The corresponding Poisson tensors P' ()2.8)1 generate elliptic and 
parabolic coordinates, tensor P' ()2.22|) gives rise to the polar coordinates and numerical 
skewsymmetric tensor generates cartesian web on the plane [2]. 

Remark 4 On the plane Q = M? there are many other solutions of the equations 1)2. 312. 4() in 
the form (E^J) 

where entries of the matrices A,B,C depend on coordinates q. However if Hi = T + V is a 
natural Hamiltonian (|1.2j) there are only four nontrivial solutions of the extended system of 
equations ()2. 312.412.12)) associated with the four known coordinate systems. 



3.2 Linear Poisson tensors and systems with cubic integrals of 
motion 

In this section we consider the following ansatz for the Poisson tensor 



P' 



( h{q) Cipi + C2P2 C3P3 + C4P4\ 

* C5P1 + CqP2 CjPi + C8P2 

* * f{q) 

* * * 



(3.6) 



Here Ck are undetermined coefficients and /, h are some functions on coordinates q. 

Substituting P' (13.6)1 into the Jacobi equations ()2.3|2.4)) we can extract equations for 
the function h. Solving these equations it is easy to prove that either function h{q) is a 
constant h{q) = const, or it is a linear function of coordinates h{q) = aqi + bq2. 

Proposition 1 // the Poisson tensor P' has the form \3. 6)) and integrals of motion 

Hi=pI+pI + Viqi,q2), H2 = H2ipi,P2,qi,q2), 

satisfy to the Lenard-Magri recurrence relations ^2.14]) then the extended system of equa- 
tions h2.S\2. 4^2.1^) with control matrix F h2.1^) has only one nontrivial solution up to 
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canonical transformations 



P' 



(0 


a 


Pi 





\ 


* 












* 








M- 

oe a 


91 


V 


* 


* 





/ 



Hi=pi+pi + 2abe'-^ . (3.7) 



The corresponding second integral of motion H2 = {pi +p2)Hi — \{pi +^2)'^ is formally 
a third order polynomial in the momenta. 

This integrable system coincides with open Toda lattice of A2 type. 

Remark 5 From the Lenard-Magri relations 1)2. 1412. 12() with control matrix F 1)2. it follows 
that 

N*dHi = dH2, N*dH2 = a2dHi + aidH2, 

N*'dH, = dHf\ N*'dHf^ = af^dH, + af^ dHf\ (^.S) 



where dH!^^ = a2dHi + a2dH2 and cr|^'* are symmetric polynomials associated with N*"^. Similar 
relations (jSIHI) may be constructed only for the A^*^ . 

In our case (|3.7() H^^ = {piP2 — ab exp((g2 — Qij/o) )^ — ^-f^i is a fourth order polynomial 
in the momenta. Below we will use deformations of the operator A^^ in order to construct open 
Toda lattices of BC2 and P2 type. To get the Toda lattice of Q2 type we can use deformations 
of the operator A^^. 

Remark 6 The Poisson tensor P' (|c{.7j) admits a natural multi-dimensional generalization 

^ dpi+i dpi ^ * dqi dp, 2 ^ dqj dq, 

which was considered in [^. The corresponding open Toda lattice of An type are one of the 
well-studied bi-hamiltonian systems. 

Throughout the rest of the section we will use the matrix F ()2.19|) associated with 
the homogeneous Stackel matrix S ()2.18|) of the Veronese type. 

Proposition 2 // the Poisson tensor P' has the form \3. 6|) and 

H,=pI+pI + V{q^,q2), H2 = //2(pi,P2,gi,g2), 

then equations \2.^2.4f^.l^) with matrix F \2.19^) has only four nontrivial solutions up 
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to canonical transformations 



P' 



P' 



P' 



/O qi -|p2 




hq\'^ 
/ 



P' 



/O -2 p2 3pi\ 

* pi P2 

* * 

* * -6e^i J 

\ 

* 

/O gi \ 

* pi -2p2 

* * 6^1 

* * / 

/O 51 4pi\ 

* pi 2p2 
6 



* * * 



- 

9i 



0/ 



V 



V 



V 



V 



--+— j +C'2e^ 



91 V^g2 " 
'2 2 , 



2/3 



2 + Ciql + C2?2 - Csg? 



Here Cj are arbitrary constants, which appear by integration of differential equations. 



The corresponding integrals of motion are the third order polynomials in the momenta p, 
which may be found in |iOj . The first solution is a periodic Toda lattice of ^3 type in the 
center of mass system ^U]. The second solution is a Holt system |1U1 ITT] . Two remaining 
solutions may be considered as integrable systems of Calogero type. 

In fact we proved that all these systems are the uniform Stackel systems with the 
homogeneous Stackel matrix S ()2.18|) . It is easy to check that the corresponding separated 
equations give rise to the hyperelhptic curves. 



Remark 7 There are many other solutions of the equations (|2.3I2.4|) if we suppose that coef- 
ficients Cfc in P' ()3.6|) are some functions of coordinates q. However if Hi = T + y is a natural 
Hamiltonian (|1.2)) there are only four nontrivial solutions of the extended system of equations 

(I2.;^I2.4I2.12I) . 

We fix the form of the desired Poisson tensor P' ()3.6|) and, therefore, in the equations 
(I2.12|l we could use the natural Hamilton function defined up to canonical transforms 



H = T + V{q), T bij{q)p,pj + ^bk {q)pk 



(3.10) 



k=l 



Here bij{q) and bk{q) are functions of coordinates q, such that kinetic energy T may 
be reduced to the canonical form ()1.2|) by point transformations q ^ q' = Aq + b and 
canonical shifts of the momenta Pi ^ Pi + bi{qi). 

As an example we consider natural Hamiltonians with pseudo-Euclidean metric jSj and 
reproduce one Drach system which is not superintegrable Stackel system with quadratic 
integrals of motion fl^. Only for this Drach system the separated variables was unknown. 
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Proposition 3 // the Poisson tensor P' has the form / Ic/. 6|) and Hi = piP2 + V{qi,q2) 



then one of the solutions of the equations 
\ 



|) reads 



P' 



/O qi 



pi + ap2 







b 

/ 



rr 2agi - q2 Ci 
Hi = piP2 - 2b — h 



C2 ( 52 - y^i 



The second integral of motion H2 is a third order polynomial in momenta. The corre- 
sponding separated variables are the special Darboux-Nijenhuis variables 



Ai,: 



P2 



± 



pI + 16b Jlfi 



4 4 

and the separated equations give rise to hyperelliptic curve. 



Remark 8 In this section we identify the first integral Hi from the Nijenhuis chain (|2.12j) with 
the natural Hamiltonian H = T + V . Of course, we can choose any integral H^ from the chain 
as the natural Hamilton function. 



3.3 The Poisson tensors second order in the momenta 

In this section we consider the following ansatz for the Poisson tensor 



P' 



( hi{q)pi + h2{q)p2 E 4jPiPj + gi(^) E 4jP^Pj + §2(9) \ 

* E4PiPj + g3(9) E4w + S4(g) 



V 



fi{q)pi + f2{q)p2 
* / 



(3.11) 

Here fi,gi,hi are some functions on coordinates q and ^ are undetermined numerical 
coefficients. 

As above there are solutions with potentials V separable in physical coordinates qi^2, 



dq\dq2 



0, and solutions with formally fourth order integrals of motion for which H2 



consists of Hi and another independent quadratic polynomial H2, for instance H2 = 
aHf + H2 or H2 = H1H2. In contrast with the previous section we will not present such 
solutions below. 

Example 7 By using tensor P' ()3.11|) we can obtain integrable systems with quadratic 
integrals of motion. For instance one of the solutions of the equations (|2.3l2.4l2.12j) with 
the control matrix F (j2.19p is equal to 



Hi 



P?+P2 + y(^1+^2 



Ci — C2 



C1 + C2 



2{qi + q2Y 2{qi - ^2)' 



and 



P' 



( 

* 



P\(l2 



P2q\ 


* 



-P? + gi(?) 

-P2P1 + g3(g) 




-P2P\ + g2(g) \ 
-pl^u{q) 
Pifiiq) +P2f2{q) 



(3.12) 
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where 



qi{Ciqf+3C2qlq2+3Ciqiq1+C2q^) rr - n <t ( f - r,-^. 



(91-92)^(91+92)^ 



92^1 gi(g) , /2 = gi gi(g) 



9i(g29j+3Ci9i92+3C29i9i+Ci9|) „ — „ r,-^ ^ i r,\ f — ^-^r 



(91-92)^(91+92)^ 



Q2qi g2{q) , fi = qi g2(g) 



The corresponding second integral of motion H2 = H2 — \ Hl is a fourth order polynomial 
in momenta, which may be reduced to a quadratic integral. 

The associated with P' ()3.12|) solution of the equations ()2.3I2.4I2.12|) with another 
control matrix F fl2.13|) has a form 

rr 2,2 Ci{ql + ql) 2C2gig2 



(qi + q2Y{,qi - q2Y (gi + q2Y{,qi - q2Y ' 
As above second integral of motion is reduced to a quadratic polynomial in the momenta. 

Integrals of motion fourth order in the momenta. 

Below we consider integrable systems for which second integrals of motion are the fourth 
order polynomials in the momenta, which can not be reduced to quadratic polynomials. 

Substituting P' ()3.11|) into the equations ()2.3|2.4|1 we can extract equations for the 
functions hi. Solving these equations we get that either functions hi are constant hi = 
const, or they are linear functions of coordinates hi = aiQi + 64^2 • 

At hi = const we consider only the Lenard-Magri recurrence relations. Recall that in 
this case the control matrix F in ()2.12j) has a special form ()2.13|) . 

Proposition 4 If the Poisson tensor P' has the form there are only two solutions 

of the extended system of equations ll2.,%2.4j2.1S^} with a control matrix F h2.1^) and a 
natural Hamilton function Hi = pl+ + V{q). 

The corresponding Poisson tensors are distinguished by functions fi and g, only 



P' 



( 2api + 2hp2 pi + gi{q) g2{q) \ 

* g-i{q) pl + gi{q) 

* * fiiq)pi + f2iq)p2 

y * * * J 



(3.13) 



If a ^ and b ^ then the Hamiltonian reads as 

Hi = pI + pI + {a + b)C2 exp ( — — ^ ) - (a - b)Ci exp ' ~^ 



a + b J \ a — b 
whereas functions fi and gi in are equal to 

/i = €26 -+b - Cic , gi = 26/2, g2 = bfi - a/2 

/2 = €26^^ + CiC , g3 = -g2, g4 = 2a/i . 

If a = or b = 0, then this is another nontrivial solution. At a = and 6=1 the 
corresponding Hamiltonian reads as 

Hi-p,+p, + 2Cie +2C2e 2{C2e^'i^ - C,f {C2e^q. - C^f ^^'^^^ 
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whereas functions fi and gt in |X23) o,re equal to 



■11 



2/2 



^1 + 



/2 = (C2e'?^ + Cle-'?2)e-'^^ 



1) 



(C2e292-Cl)2 (C2e292-Ci)2 



In the both cases the second integrals of motion are the fourth order polynomials in the 
momenta. 

The first solution is an open Toda lattice of P2 type. If C3 = C4 = the second 
Hamiltonian Hi flH.14|l describes open Toda lattice of type. If Ci = or C2 = the 
Hamiltonian Hi flH.14|l describes also an open Toda lattice associated with a root system 
BC2- If Ci = C2 we obtain the Hamiltonian Hi fl3.14|) of the generalized Toda lattice, 
which was constructed by Inozemtsev 13j. 

Remark 9 The quadratic matrix polynomial P' 1)3. may be considered as deformation of 
the following product of linear Poisson tensor (|ll7j) (see Remark 



P' — PlinP ^Plin + ^lin ; 



(3.15) 



where Knj^^ is a very simple linear matrix polynomial. Therefore, substituting the multi- 
dimensional linear Poisson tensor (|3.9|) into 1)3. 15() and adding the similar term one gets a 
natural multi-dimensional generalization of P' ()3.13|) . which was proposed in j4j for the open 
Toda lattices at C3 = C4 = and at Ci = (C2 = 0). 

Now let us consider the second case at hi = a^qi + biq2. 



Proposition 5 // the Poisson tensor P' has the form there are only two explicit 

nontrivial solutions of the equations \2.'J\2.4\2.1^/3^] with control matrix F \2.19\i and a 
natural Hamilton function Hi = p\+ p\ + '^{o)- 

The first solution is the Holt-type system for which 



Hi=pi+pi- 



6q. 



2/3 



Ql q^ 



and 



/o 



V 



3pig2 




3C2g2 



<ii 




pj 



P1P2 + 61^1^2 
3Ciqt^' + ^ + 2Csqf- 
-Cipiql^' - ^ 

^ J 



The second solution is the system with fourth order potential for which 



Hi=pI+pI 



C.(4,^ + 3,?.,^ + |)+C2(,?4)-|-- 



(3.16) 



(3.17) 



(3.18) 



and 



P' 



/O piq2 p? + gi(g) piP2 + g2{q) \ 

* g3(g) P? + 2p^ + g4(g) 

* * +^2/2(5) 
y* * * 



(3.19) 
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where 



C 

g2 = CiQiql, gi = -4Cig^ + C2ql 1, 



9?' 



g3 = 4g2 + 92/2, g4 = -C,{Aqt + 6gfgi + g|) + C2 ( gf + ^''^ 



r/ie second integrals of motion are the fourth order polynomials in the momenta llUf . 

According to in the second case we can make a special contraction of variables and 
obtain the Henon-Heiles system with the cubic potential 



H,=pI+pI + igi(16g2 + 3^2') + ^ . (3.20) 
"J Q'2 



The corresponding Poisson tensor P' is equal to 



P' 



( 92^1 pI + y g? P1P2 - \ql \ 

* 8g2g? + |gf pi + 2pl + 2q,ql + fqf 

* * -^2(851+52) 
* / 



(3.21) 



Remark 10 The eigenvalues of the operators N associated with the Poisson tensors ()3.19() and 
(|H.21|) are integrals of motion Ai,2 = 0. So, the special Darboux-Nijenhuis coordinates can not 
be the separated variables for the system with fourth order potential H3.18() and for the Henon- 
Heiles system (|ll'2()|) . The corresponding separated variables are generic Darboux-Nijenhuis 
coordinates {x,y) 1)2. 9|) . 

If the Poisson tensor P' has a form 1)3.111) then there are some implicit solutions of the 
equations ()2. 312. 412. 12)) with control matrix F ()2.19|) . In these cases functions fi{q),gi{q) 
and integral of motion H2{p, q) are determined by potential ^(^1, ^2), which satisfies some 
special partial differential equations of the second and of the third order. For instance, 
one of the obtained equations has the form 



Similar second order equations were considered in ^01 . The third order partial dif- 
ferential equations obtained in this method we could not find in the literature. These 
imphcit systems will be studied separately. 

As above we can substitute into the equations ()2.12p the Hamilton function Hi = T+V 
in the form ()3.1(J|) . Here we present two new integrable systems with pseudo-Euclidean 
metric. 

Proposition 6 // the Poisson tensor P' has a form VJ.ll]) there are only two explicit 
solutions of the equations i2.^2.4^2.1^} . The first solution reads as 

Hi = P1P2 + Ci(4g2 + 5agig2 + a\l) + ^2(91 + 0^92) ^ _^ ^ _^ ^ ^ ^ 

92 92 92 
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/ 3(72P1 pf-6Cig2(qi+ag2) « - 3pip2 -3Ci (gi +aq2 ) (2gi +ag2 ) - '^^ '"1/3°"^ ' + ^ \ 

«2 «2 



P' 



ISCig^ 




2/3 



9Cig2(2gi+Q:g2)+3C2g2 

3Ci(2gi-ag2)+-^+ )pi+18Ci 




192P2 



The second solution has the form 



Hi = P1P2 + ( Ci + — L ) qi + C-iq2 H t^v^ , 



4g2Pi 



Pl+'7l+0'72 
252 


* 



(6C1 +4a)pf - 2pip2 +2(Ci +a)gi +2Ci a ( (^+ ) ' „ ^i^^||±^) \ 
P?+gi+{a+4Ci)g2+4C2V92 

-CTP1+P2 

/ 



where a 



C:i-2Cf 



In the both cases the corresponding second integrals of motion are the fourth order poly- 
nomials in the momenta. Moreover, the special Darboux-Nijenhuis coordinates associated 
with are the separated variables and the separated equations give rise to hyperelliptic 
curves. 

3.4 The normalized Poisson tensors 

In this section we consider the normalized Poisson tensors 

P' = det(P')"^P', (3.22) 

where P' has the form (|3.11|) . Such tensors induce formal ujN structure on M, but their 
eigenvalues Aj will be a priori functionally dependent functions. 



Proposition 7 For the Poisson tensor P' 2^) there are only two nontrivial solutions 
of the equations i2.^2.4^2.1^} with the control matrix F \2.19jl . 



The corresponding Poisson tensors are distinguished by functions fi and gi only 



P' 



( v\(i2 -pi + gi(g) Pij52 + g2(g)\ 

* g3(g) P2 + g4(?) 

* * fiiq)pi + f2iq)P2 

y * * * J 



(3.23) 



First solution corresponds to the Henon-Heiles system for which 

Hi=p\+pI- CMlQql + Ml) - 2C2qi, 
whereas functions fi and gi in \3. 2!^) are equal to 

fi = 0, gi = SCiqiql g2 = ^ql 

/2 = 3Ci(8gf + g|) + C2, g3 = 4g2 - /2g2, g4 = -gi • 
Second solution corresponds to the system with the fourth order potential for which 

H,=pI+pI + Ci{8qt + Qqlql + qt) + C2{^1 + ql) + 
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whereas functions fi and gj in (Eil^j are equal to 

h = 2C,ql gi = -{Ci{ql + 6q!) + C2)ql = -2C,q,ql 

f2 = -4Cigi (3gi + Aqf) - 4^2^! + ^, gs = 4g2 - /aga, g4 = -gi • 

The corresponding second integrals of motion are equal to 

H2 = det(P') = ^Jpi + qifipl + (^2/2 + gz)PiV2 - 2gipi - g? + gags • 

In the both cases eigenvalues of the recursion operator N are trivial Ai 2 = ±1. In order 
to get the separated variables we have to put operator in diagonal form ()2.9|) by using 
canonical transformation z = {q,p) ^ Q = {x,y). 

Proposition 8 Let J = %. is the Jacobi matrix associated with some change of variables 
C = C(-2)- Solving overdetermined system of the partial differential equations 

JPJ^ = P, JA^J-i = diag(l, -1,1,-1) (3.24) 

with respect to ({p, q) we reproduce the separated variables for the Henon-Heiles system and 
the separated variables for the system with fourth order potential, which were constructed 
in H^. 

The first equation in ()3.24|) provides canonicity of transformation z ^ (, the second 
one provides diagonability of in ^ variables. Equations ()3.24|) were easy solved in the 
symbolic computational system Maple. 

Relations of such transformations, which put operator in diagonal form, with the 
Backlund transformations for these systems will be discussed in forthcoming publications. 

4 Appendix: a programm for the search of natural 
integrable systems 

This paper as well as [23 121] belongs mostly to computational mathematical physics and, 
therefore, in this section we present an implementation of the discussed algorithm made 
in the symbolic computational system Maple. 

At the first step we need to determine dimension of the phase space M and canonical 
variables g = (gi, . . . , g„) and p = {pi, . . . ,pn): 

> restart: with(linalg) : 

> n:=2: 

> q:=seq(q| I i,i=l . . n) : p : =seq(p I I i , i=l . . n) : 

> var:=q,p; 

Now we can introduce canonical Poisson tensor P ()2.6|) 

> ed : =array (identity , 1.. n,l.. n) : z : =array (sparse , 1 . . n,l.. n) : 

> PI :=blockmatrix(2,2, [z,ed,-ed,z] ) ; 

and a linear ansatz for the second Poisson tensor P' ()3.6p at h{q) = aqi 

> P2 := array (antisymmetric , 1 . . 2*n,l.. 2*n) : 

> P2[l,2] :=a*ql: P2 [1 , 3] : =cl*pl+c2*p2 : P2 [1 ,4] : =c3*pl+c4*p2 : 

> P2[2,3] :=c5*pl+c6*p2: P2[2,4] :=c7*pl+c8*p2: P2 [3 ,4] : =f 1 (q) : 

The Poisson brackets associated with the Poisson tensor P' and with the Poisson pencil 
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P + \P' are defined by the standard rule 



> P3:=evalm(Pl+lambda*P2) ; 

> PB2 := proc (f , g) options operator, arrow: 

> add(add(P2[i,j]*diff (f ,var[i])*diff (g,var[j]) ,i=l. . 2*ii),j = l.. 2*ii) 

> end: 

> PBS := proc (f , g) options operator, arrow: 

> add(add(P3[i,j]*diff (f ,var[i])*diff (g,var[j]) ,i=l. . 2*n),j=l.. 2*n) 

> end: 



The corresponding Jacobi identities ()2.2|) look like 

> f:=F(var); g:=G(var): h:=H(var): 

> Y2:=simplify(PB2(PB2(f ,g) ,h)+PB2(PB2(g,h) ,f )+PB2(PB2(h,f ) ,g)) : 

> Y3:=simplify(PB3(PB3(f ,g) ,h)+PB3(PB3(g,h) ,f )+PB3(PB3(h,f ) ,g)) : 

It allows US to build up the system of equations ()2.3I2.4|) 

> ListEq:=NULL: 

> for i to 2*n do 

> for j to 2*n do 



> for k to 2*n do 

> Eq2 : =coef f (coef f (coef f ( 

> Y2,diff (f ,var[i]) ,1) ,dif f (h, var [j] ) , 1) ,dif f (g, var [k] ) , 1) : 

> Eq3 : =coeff (coef f (coef f ( 

> Y3,diff (f ,var[i]) ,1) ,dif f (h, var [j] ) , 1) ,dif f (g, var [k] ) , 1) : 

> ListEq:=ListEq,coeffs(collect(Eq2,{p}, distributed) ,{p}-) , 

> coef fs(collect(Eq3,{p, lambda}, distributed) ,{p, lambda}) ; 



> end do : end do : end do : 

On the next step we can construct the operator A^* and its minimal polynomial A(A) 

> N:=evalm(inverse(Pl)&*P2) ; 

> ed2 :=array (identity, 1. . 2*n,l.. 2*n) : 

> Delta: =collect (simplify (sqrt (f actor ( 

> det (N-lambda*ed2) ) ) , symbolic) , lambda, factor) ; 

> if coef f (Delta, lambda, n) <0 then Delta:=-Delta: end if: 

The coefficients of the minimal polynomial A(A) are entries of the control matrix F (j2.19|) 

> Sigma : =vector (n) : 

> for k to n do sigma[k] :=-coeff (Delta, lambda, n-k) ; end do: 

> F : =delcols (augment (sigma, ed) ,n+l . . n+1) : 

Now we introduce natural Hamilton function H = Hi and unknown integrals of motion 
H2, ■ ■ ■ , Hn and calculate the corresponding differentials dH^ 

> Hl:=add(p| |k~2,k=l. . n)+V(q); 

> for k from 2 to n do Hi |k:=h| Ik(var) ; dH| |k:=vector(2*n) : end do: 

> for k to n do 

> for i to 2*n do dH| |k[i] :=simplify(diff (Hi Ik, var [i] ), symbolic) : 

> end do : end do : 

On the final step after building of the equations (j2.12|) 

> ListEqH:=NULL: 

> for i to n do 

> Z| |i:=simplify(evalm(N&*dH| |i-add(F[i,j]*dH| I j,j=l. . n) ) ) : 

> ListEqH:=ListEqH,seq(Z| |i[k] ,k=l. . 2*n) : 

> end do : 

we solve the complete overdetermined system of algebro-differential equations 

> Ans:=pdsolve({ListEq,ListEqH},{f l(q) ,V(q) ,seq(H| |k,k=2. . n)}. 
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> 



parameters={seq(c I I j , j=l . . 8) }) ; 



This programm allows us to construct three integrable systems listed in the Proposition 
121 in a few second on the standard personal computer. A little large programm allows us 
to construct all the systems listed in this paper in a few minutes. 

Remark 11 In order to cut the computer time and to exclude the trivial solutions we have 
to add some simple equations to our system. For instance we added the following inequalities 

5 Conclusion 

Using some recent results on the separation of variables for bi-hamiltonian manifolds 
we propose a new method for construction on the natural integrable systems. In this 
method the desired integrals of motion are solutions of the overdetermined system of 
algebro-differential equations, which arise from the Jacobi identities for the Poisson pencil 
(|2.3l2.4j) and from invariance of the corresponding bi-Lagrangian distribution with respect 
to recurrence operator (j2.12j) . The algorithm may be easy realized in any modern 
symbolic computational system. One of the possible implementations may be found in 
the Appendix. So the search for natural integrable systems is reduced to simple computer 
calculations. 

The main advantage of the proposed method is that we construct integrals of motion of 
a natural integrable system, the corresponding Poisson pencil and a part of the separated 
variables simultaneously. 

The main disadvantage is that we have to use undefined polynomial anzats for the 
second Poisson tensor P' and to choose some control matrix F. So, we construct special 
families of the integrable systems associated with the a given P' and F. 

Of course, by using very simple substitutions (j3.6l3.11|) we reproduce a majority of 
the known natural integrable systems on the plane with cubic and quartic integrals of 
motion. Nevertheless we did not get all the known integrable systems ^U]. Moreover we 
do not know how to construct similar polynomials P' on the sphere or on the ellipsoid. 

The second disadvantage is an exponential growth of the computer time with the 
growth of n = rankP' and with the growth of the higher power of the matrix polynomial 
P'. 

The author want to thank I.V. Komarov and M. Pedroni for valuable discussions. 

References 

[1] C. Bartocci, G. Falqui, M. Pedroni, A geometric approach to the separability of the 
Neumann-Rosochatius system, Differential Geom. AppL, v. 21(3), p. 349-360, 2004. 

[2] S. Benenti, Orthogonal separable dynamical systems. In Proceedings of the 5th In- 
ternational Conference on Differential Geometry and Its Applications, Silesian Uni- 
versity at Opava, August 24-28, 1992, O. Kowalski and D. Krupka Eds. Differential 
Geometry and Its Applications, v.l, p. 163-184, 1993. 

[3] R. Brouzet, About the existance of recursion operator for completely integrable 
Hamiltonoan systems near a Liouville torus, J. Math. Phys., v. 34, p. 1309-1313, 1993. 



20 



[4] P.A. Damianou, Multiple Hamiltonian structures for Bogoyavlensky-Toda lattices, 
Reviews Math. Phys., v.l6 p.175-241, 2004. 

[5] J. Drach, Sur I'lntegration logique des equations de la Dynamique a deux variables, 
Comptes Rendus, Paris, v.200, p.22-26, 1935. 

[6] G. Falqui, F. Magri, M. Pedroni, Bihamiltonian geometry and separation of variables 
for Toda lattices, J. Nonlinear Math. Phys., v.8, p. 118-127, suppl. (Proceedings of 
NEEDS99), 2001. 

[7] G. Falqui, M. Pedroni, Separation of variables for bi- Hamiltonian systems, Math. 
Phys. Anal. Geom., v.6, p.139-179, 2003. 

[8] A. Fokas, B. Fuchssteiner, Symplectic structures, Bdcklund transformations and 
hereditary symmetries, Physica D, v.4, p.47-66, 1981. 

[9] R.L. Fernandes, On the master symmetries and hi- Hamiltonian structure of the Toda 
lattice, J. Phys. A, v.26, p.3797-3803, 1993. 

[10] J. Hietarinta, Direct method for the search of the second invariant, Phys. Rep., v. 147, 
p.87-188, 1987. 

[11] C.R. Holt, Construstion of new integrable Hamiltonians in two degrees of freedom, 
J. Math. Phys., v.23, p.1037-1046, 1982. 

[12] A. Ibort, F. Magri, G. Marmo, Bihamiltonian structures and Stdckel separability, J. 
Geometry and Physics, v.33, p.210-228, 2000. 

[13] V.I. Inozemtsev, The finite Toda lattices, Commun. Math. Phys., v.l21, p. 629-643, 
1989. 

[14] Y. Kosmann-Schwarzbach, F. Magri, Poisson-Nijenhuis structures, Ann. Inst. 
Poincare (Phys. Theor.), v.53, p.35-81, 1990. 

[15] T. Levi-Civita, Integrazione delle equazione di Hamilton- J acobi per separazione di 
variabili. Math. Ann., v.24, p.383-397, 1904. 

[16] F. Magri Geometry and soliton equation. In: La Mecanique de Lagrange et son 
heritage, Atti Acc. Sci. Torino Suppl, v. 124, p. 181-209, 1990. 

[17] F. Magri, Eight lectures on Integrable Systems, Lecture Notes in Physics, Springer 
Verlag, Berlin-Heidelberg, v.495, p.256-296, 1997. 

[18] C. Morosi, G. Tondo, Quasi-Bi-Hamiltonian systems and separability., J. Phys. A., 
V.30, p. 2799-2806, 1997. 

[19] C. Morosi, G. Tondo, The quasi- bi- Hamiltonian formulation of the Lagrange top. J. 
Phys. A., V.35, p.1741-1750, 2002. 

[20] S. Rauch-Wojciechowski, A.V. Tsiganov, Quasi-point separation of variables for 
Henon-Heiles system and system with quartic potential, J. Phys. A., v. 29, p. 7769, 
1996. 

[21] A.V. Tsiganov, The Stackel systems and algebraic curves, J. Math. Phys. , v.40, 
p.279-298, 1999. 



21 



[22] A.V. Tsiganov, On the Drach superintegmble systems, J. Phys. A., v. 33, p. 7407, 2000. 

[23] A.V. Tsiganov, Yu.A. Grigotyev, Symbolic software for separation of variables in the 
Hamilton- Jacobi equation for the L-systems, Regular and Chaotic Dynamics, v. 10(4) 
p.413-422, 2005. 

[24] A.V. Tsiganov, Yu.A. Grigotyev, Computing of the separated variables for the 
Hamilton- Jacobi equation on a computer, Nonlinear mechanics, v. 1(2), p. 163-179, 
2005. 

[25] M. Wojciechowska, S. Wojciechowski, New integrable potentials of two degrees of 
freedom, Phys. Lett. A., v.l05, p.11-14, 1984. 



22 



